Effects of Anisotropic Charge on Transverse Optical Phonons in NiO 
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Phonon dispersion of detwinned NiO is measured using inelastic x-ray scattering. It is found 
that, near the zone center, the energy of the transverse optical phonon mode polarized parallel to 
the antiferromagnetic order is ~1 meV lower than that of the mode polarized perpendicular to the 
order, at room temperature. This is explained via anisotropic polarization of the Ni and O atoms, as 
confirmed using a Berry's phase approach with first-principles calculations. Our explanation avoids 
an apparent contradiction in previous discussions focusing on Heisenberg interaction. 

PACS numbers: 75.50.Ee, 75.30.Gw, 77.84.Bw, 63.20.dd 



Transition-metal mono-oxides are fundamental materi- 
als for studying properties of magnetic and strongly cor- 
related systems. Among them, NiO and MnO are antifer- 
romagnetic (AFM) insulators with similar crystal struc- 
ture and physical properties. The AFM order and lat- 
tice contraction occur along the [111] direction below the 
Neel temperature (Tn=523 K for NiO). The order con- 
sists of ferromagnetic (111) planes. Above Tn, both have 
a rock-salt structure. NiO and MnO are relatively sim- 
ple, because they have non-degenerate electronic ground 
states, and are free from Jahn- Teller effects, unlike FeO 
and CoOii 

The magnetism of these materials is often discussed in 
terms of a superexchange mechanism. MnO has antifer- 
romagnetic interactions both in the nearest neighbor ( Ji) 
and next nearest neighbor (J2) exchanges, and these in- 
teractions reproduce the experimental results.^'"^ On the 
other hand, Ji for NiO still remains uncertain, perhaps 
because it is much smaller than the antiferromagnetic su- 
perexchange J2 interaction, which is responsible for the 
AFM order. For example, local spin density approxi- 
mation (LSDA)-fJ7 calculations give ferromagnetic Ji^ 
while some calculations such as GW show antiferromag- 
netic Ji.** Experimentally, a measurement of spin- wave 
dispersion indicates ferromagnetic Ji,— while one of mag- 
netic susceptibility suggests antiferromagnetic Ji.=£ 

Another issue is that LSDA-I-J7 calculation results are 
not easily reconciled with the exchange interaction pic- 
ture of NiO. According to Refsi^ii, the lattice distortion 
is dominated by Ji with no contribution from J2, and 
antiferromagnetic Ji causes the contraction in [111], if 
|Ji| decreases with increasing the distance between the 
nearest-neighbor Ni atoms. Based on this discussion, the 
calculated ferromagnetic Ji^ is not consistent with the 
calculated^ (and observed) contraction in NiO, suggest- 
ing an additional ingredient is needed to understand the 
calculation results. 

The energy of transverse optical (TO) phonons can 
be used as a direct probe of microscopic coupling^; TO 
modes that would be degenerate in a cubic rock-salt 
structure, are split at the zone center under the AFM 
order, with the energy of the mode polarized along the 
order (-£'™) different than that of the mode polarized 



in the plane perpendicular to the order (£'™). As dis- 
cussed in Ref.'^ the sign of this difference (i?™— iJj*-*) 
can be linked to the sign of Ji. The picture is gener- 
ally consistent for MnO, where the observed shif t^°'^^ is 
in agreement with theory^^i^ however question remains 
about NiO, where preliminary measurementa^^ disagree 
with the hST)k+U calculations 

We have made careful measurements of the TO phonon 
dispersion in detwinned NiO at room temperature (RT), 
and find i?™— ijj^^— 1.0 meV in the vicinity of the 
zone center, in reasonable agreement with the LSDA-t-J7 
calculations.^'® In order to understand this, we intro- 
duce the anisotropic polarization derived from a Berry's 
phase.— This anisotropy yields an anisotropic effective 
Born charge, which was introduced into AFM materials 
by Ref.^, after tremendous success in ferroelectric mate- 
rials. This spin-dependent term dominates the splitting 
and provides the missing ingredient to explain the ob- 
served lattice distortion and phonon splitting simultane- 
ously. 

NiO samples with (111) surface were detwinned in ac- 
cordance with Ref.— , giving the AFM order and lattice 
contraction along the [111] direction (in this report, in- 
dices of cubic symmetry are used). The single structural 
domain (T-domain) was confirmed by the (333) Bragg 
peak (Fig. 1(a)); a peak from another domain with a dif- 
ferent direction for the contraction, if any, would appear 
at lower angle (the bar in Fig. 1(a)). NiO may also have 
a different type of domains, originating from the differ- 
ent spin-axis orientation (5-domain). However, Ref.— 
shows that a field of >2 kOe along [110] may be used to 
align the spin axes along [112]. When the field reaches 5 
kOe, a single S'-domain is obtained,— but the T-domain 
wall starts to move.--" Considering these facts, magnetic 
field of 4.5 kOe was applied parallel to [110] during the 
measurements. 

Phonon measurements using inelastic x-ray scattering 
(IXS) were performed at BL35XU of SPring-S^^ with a 
total energy resolution of ^1.6 meV and a momentum 
resolution of ^0.08 A~^. The incident energy is 21.75 
keV. The measurements were achieved at RT (^300 K 
(<T/v)) to prevent introducing a new domain caused by 
tiny stresses or thermal gradients J^ii^ We also performed 
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FIG. 1. (Color Online) (a) The 9 scan around the (333) Bragg 
peak. The bar indicates the expected peak position of {333}. 
(b) IXS spectra at \q\ — {0.1, 0.1, 0.1). LA denotes longitu- 
dinal acoustic phonon. TAs (LAs) indicates the Stokes line 
of TA (LA), while TAa (LAa) does the anti-Stokes line, (c) 
The enlarged view of (b) at around the TO phonons. Dotted 
lines are results of curve fitting, and the bars show the peak 
positions. 



similar measurements under 1.7 kOe at some phonon vec- 
tors to study field induced effects, such as magnetoelec- 
tric effects. The results show the same phonon peak po- 
sitions as at 4.5 kOe, within an experimental accuracy of 
±0.3 meV. We confirmed the zero energy using the mid- 
point of the Stokes and anti-Stokes lines of the transverse 
acoustic (TA) phonons. 

Figure l(b, c) shows the IXS spectra at a phonon vec- 
tor of |(7| = (0.1, 0.1, 0.1), where q—Q—r; Q is the total 
momentum transfer, and r is the nearest Bragg point. A 
shift is observed, depending on the choice of momentum 
transfer. For the top spectrum, where q is parallel to 
[111], the polarization vector of the observed transverse 
mode must be in the plane perpendicular to the [111] 
ordering direction. For the others, the situation is less 
well defined, and we will discuss it again in the context 
of specific models. However, naively, one would expect 
the polarization vectors to be within and perpendicular 
to the plane defined by the [111] axis and q. In this case, 
at the chosen total momentum vectors, Q, the mode po- 
larized in the plane would tend to dominate the lower 
three spectra. 

Simulations were carried out using shell model 
calculations,^^ which describe phonon dispersion, with 
only a few parameters, very well especially for NiO and 
MnOjiiii^ The OPENPHONON source code^^ was used for 
the calculations, and most of the parameters were taken 
from Rcf.^''.^'^ The shell model has several types of pa- 
rameters including the atomic position. Coulomb inter- 
action (charge), including the core-shell interaction, and 
short-range Born-von Karman force constants between 
pairs of atoms. As seen in Fig. 2(a), the calculations as- 
suming the cubic structure reproduce the present exper- 
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FIG. 2. (Color Online) (a) Phonon dispersion of NiO in the 
|g| = [qqq] direction, which is determined theoretically (lines), 
and experimentally (the dots) . (b) Simulated TO phonon dis- 
tribution in the vicinity of |Q| = (333) (solid line). Dotted line 
shows the spectrum broadened by a 0.4 meV resolution func- 
tion, (c) The band structure along [qqq] and [qqq] obtained 
by LSDA-ft^. Solid lines are the valence bands which are used 
in calculating the Berry's phase of polarization. 



imental phonon dispersion very well. This agreement in 
dispersion lends credence to the idea that our approach, 
below, of modifying successively, the various components 
of the shell model, can be used to understand the mate- 
rial's behavior. 

First, just the effect of the non-cubic lattice was con- 
sidered, holding all the other parameters fixed. The ex- 
perimental lattice distortion of ~90.06° at RTi^ gives TO 
phonon splitting of -0.06 meV at |gl=(0.1, 0.1, 0.1). This 
value is much smaller than the observed in Fig. 1(c), sug- 
gesting that only the lattice distortion can not describe 
the observed TO peak shift. 

Next we introduce additional Born-von Karman force 
constants to simulate the anisotropy of the Heisenberg 
exchange interaction. Effectively, the shell model in- 
cludes J2 in the longitudinal force constant between near- 
est neighbor Ni-0 atoms. As for Ji, additional force con- 
stants along (Oil) (between parallel spin, T) and (Oil) 
(between antiparallel spin, —T) were included (as the 
cubic structure assumed, Ji deviation due to the lattice 
distortion^ (>^i^) is not considered). The parameter T 
was varied to fit the experimental results. 

Using these additional force constants, the polarization 
vectors for the peaks in Fig. 1(c) can be determined as 
eigenvectors of the dynamical matrix. The spectra at 
q=(0.1, 0.1, -0.1) and (0.1, -0.1, 0.1) with r=(333) are 
assigned to the modes polarized close to [112] (T0[ii2]) 
and [121]_(TO[i2i]) respectively. For g=(0.1, -0.1, 0.1) 
with T=(151), the spectrum is dominated by the T0ji2i] 
mode, while the spectra at q—{0.1, 0.1, 0.1) with t=(151) 
is assigned to TOj_. In short, it is estimated from the 
experiment that the energy of T0(ii2) {Ej^2)) ~0.9 
meV lower than EJ^. 

The intensity around the zone center is simulated us- 
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ing the parameters above. Figure 2(b) shows the calcu- 
lated TO lineshape (solid line) near |(5| = (333) consists 
of two parts; a delta-function-like TO_l peak and distri- 
bution from TO II to TOj^ resulting from the direction 
dependence of q^O. From Fig. 2(b), difference between 
i?™ and EJ'-' is estimated as '^1.0 meV. The value at 
RT agrees reasonably with the results of the LSDA+U 
calculations (i?™— i?™=— 1.8 meV),"^'® which generally 
assume r=0. Furthermore, when this spectrum is con- 
volved with a 0.4 meV resolution function (dotted line), 
it is clearly seen that TO^ has stronger intensity than 
other components. The stronger TOj^ intensity agrees 
with the observation in MnOJ^iii 

Both the phonon splitting and lattice contraction for 
NiO are correctly predicted by the LSDA+U calcula- 
tions. However, conceptually LSDA-|-[/ does not fit the 
superexchange mechanism well; the calculated ferromag- 
netic Ji and lattice contraction require increasing |Ji| 
with increasing r in the exchange framework (r indi- 
cates that the distance between the nearest-neighbor Ni 
atoms), since the lattice distortion is proportional to 
—dJi/dr^^ Moreover, this contradicts the relation, sug- 
gested in Ref.— , between the sign of £'™— and that 
of Ji. Thus some additional ingredient is required, which 
is included in the calculations but not in the exchange 
framework, in order to explain the experimental observa- 
tion. 

The last perturbation we consider is an anisotropic 
Coulomb interaction. The effective Born charge tensor, 
Z*, was estimated using LSDA-I-J7 as implemented in the 
ABINIT package)^ U =4 eV was chosen to fit the exper- 
imental valence band,— mainly because the calcula- 
tions for the Berry's phase of polarization only require 
the occupied states (solid lines in Fig. 2(c)). The calcu- 
lated charge tensor of Ni is 

/ z* Az* Az* \ 
Z^i = Az* z* Az* (1) 
y Az* Az* z* J 

where z*=2.26 and Az*==0.03. For oxygen (Z^), 
z*=— 2.26 and Az*=— 0.03 are obtained independently, 
satisfying the charge neutrality. This gives an anisotropic 
parameter, Az*/z*=1.3x 10~^. The principal directions 
and values of the tensor show that each atom has an 
anisotropic charge of z*-t-2Az*parallel to the AFM or- 
der, [111], (Z|p andz*-Az* perpendicular to [111] (Z^). 
The anisotropy results from the asymmetric band dis- 
persion; the AFM order brings about band folding in 
the [qqq] direction (Fig. 2(c)), while no folding oc- 
curs in the (qqq) directions, leading to the effective 
charge anisotropy. Hence this charge anisotropy is spin- 
dependent, extending, somewhat, the limits of Ander- 
son's original discussioni^i 

This anisotropic change may be introduced into the 
shell model in different ways. In particular, the detailed 
formulation of the shell model introduces matrices for 
the total charges Z, shell charges Y and core-shell force 



constants k. We could introduce anisotropy into any of 
these, Z, Y, k, and here we consider explicitly Z and Y. 

The free-ion polarizability is given by Y'^e'^/k so that 
we first consider introducing the anisotropy in Y. As- 
suming the off-diagonal components, Yij (iy^j), for both 
O and Ni have a common value of Yij /Yii^2. 8 xlO^"^, we 
can reproduce the experimental results in Fig 1(c) and 
all assignments mentioned above. Based on Ref.^*', which 
discusses the shell model in detail for a rock-salt struc- 
ture, we then estimate the effective charge anisotropy as 
Az*/z*=5.6xl0"2 (z*=±2.21). 

Alternatively, we consider introductions of anisotropy 
directly in the static charge, Z, of the shell model. Tak- 
ing Zij / Zii=3xl0~^ is sufficient to reproduce the results. 
This value corresponds to Az*/z*=8xl0-3 (z*=±2.19). 
This agrees better with the LSDA-|-[/ calculations than 
the value from the anisotropic Y case above. However, it 
is also possible the larger Az*/z* value is caused by lack 
of consideration of higher order moments than dipole in 
the shell model.— 

It is interesting to consider if this anisotropy might be 
related to the static contraction along the [111] direc- 
tion. According to Ref.—, Born effective charges have a 
correlation to the covalent bonding. This suggests the 
lattice distortion can be explained using the anisotropic 
Born charge, as a distance between two atoms at equi- 
librium may depend on the magnitude of the charges. 
Noting that the principle directions of Z* are orthogo- 
nal, the lattice distortion from cubic to rhombohedral 
can be estimated using Ref.—. In the anisotropic Z case, 
Az*/z*=8xl0~^ gives the lattice distortion of 90.15°, 
in fair agreement with the measured 90.06° at room 
temperature.^^. If we take Az*/z*=5.6x 10^^, the dis- 
tortion in the case of anisotropic Y is estimated, giving 
a large value, 91.11°. 

The measured momentum dependence of the splitting 
is in agreement with both of these (anisotropic Y or 
anisotropic Z) models as can be seen in Fig. 3 (a,b). One 
notes that in all models, the phonon splitting is reduced 
towards the zone boundary, and the magnitude at the 
zone boundary depends on the detailed model. However, 
our resolution does not allow us to choose between the 
models. Meanwhile, Figure 3(c) depicts the TO distribu- 
tion near |Q| = (333) for anisotropic Y and Z, showing the 
similar distribution to Fig. 2(b) for the anisotropic force 
constant. It should be noted that, based on Ref.^^, it 
is confirmed that these anisotropic charge models satisfy 
the generalized Lyddane- Sachs- Teller relatio n^'^i^^ both 
along and perpendicular to the [111] direction. 

In conclusion, we have two different methods of in- 
cluding the anisotropy of the material that fit our ex- 
perimental results for the phonon splitting: (1) ad hoc 
addition of force constants to mimic a possible Heisen- 
berg interaction and (2) inclusion of an anisotropy in the 
electrostatic terms of the shell model. The latter, which, 
qualitatively seems to match the anisotropy of the Born 
effective charge tensor calculated using the LSDA+f7 
method, seems a more reasonable choice. Furthermore, 
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FIG. 3. (Color Online) (a) TO phonon dispersion in the [111] 
and [111] directions. Dotted lines are results of curve fitting, 
(b) Comparison of TOx and T0[ii2] phonon modes between 
experimental (marks) and theoretical (lines) results. The cal- 
culations are based on the shell model with (solid and dotted 
lines)/without (dashed lines) effective charge anisotropy. (c) 
Simulated TO phonon distribution for the effective charge 
anisotropy in the vicinity of 1Q| = (333). 



the latter description provides an additional ingredient 
that allows one to reconcile the observed lattice contrac- 
tion with the ferromagnetic Ji interactions calculated us- 
ing LSDA-Ht/^'^ 



These discussions suggest the anisotropic polariza- 
tion should be considered in studying AFM compounds. 
Both the Berry's phase calculations of polarization and 
Heisenberg exchange interaction have spin-dependent 
anisotropy. However, the exchange interaction applies 
to only magnetic atoms, while the anisotropic effec- 
tive charge applies to both magnetic and non-magnetic 
atoms, as required by charge neutrality. In addition, in 
order to estimate the exchange interaction precisely, a 
deviation of the Wannier functions from the cubic sym- 
metry should be considered, even if the lattice is assumed 
as cubic. This may be relevant to other AFM materials; 
for example, high-Tc cuprates, or iron pnictides. 

In summary, we observed TO phonon shifts for de- 
twinned NiO at room temperature, and found that the 
energy of the TO phonon with atomic motions in the di- 
rection of the [111] magnetic ordering is ^1.0 meV lower 
than the TO mode with motions perpendicular to the 
ordering direction, near the zone center. We suggest this 
observation is most simply explained by considering the 
anisotropic polarization of the medium. While ultimately 
related to a spin dependent interaction, at a conceptual 
level this is an additional ingredient that should be con- 
sidered, to understand the material behavior, distinct 
from the exchange interaction. 

The synchrotron radiation experiments were per- 
formed at the BL35XU in the SPring-8 with the approval 
of the Japan Synchrotron Radiation Research Institute 
(Proposal No. 2008A1584). 
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